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Abstract. - The Csiszar f-divergence, which is a class of information distances, is known to offer 
a useful tool for analysing the classical counterpart of the cloning operations that are quantum 
mechanically impossible for the factorized and marginality classical probability distributions under 
Liouville dynamics. We show that a class of information distances that does not belong to this 
divergence class also allows for the formulation of a classical analogue of the quantum no-cloning 
theorem. We address a family of nonlinear Liouville-like equations, and generic distances, to obtain 
constraints on the corresponding functional forms, associated with the formulation of classical 
analogue of the no-cloning principle. 



Introduction. — The search for classical analogues of the processes that are quan- 
tum mechanically impossible allows one to identify those aspects that are purely quantum 
mechanical. Along this line, a classical analogue of the quantum entanglement has been 
discussed, namely, secret classical correlations [1]. Moreover, it has been reported that prin- 
ciples of quantum no-cloning [2, 3] and no-deleting operations [4] under unitary dynamics 
possess classical analogues [5]. The existence of parallelism between classical and quantum 
copying processes has been demonstrated via contradictions between the behaviors of rela- 
tive entropies. That is, the relative entropy calculated from statistical ensemble distributions 
changes upon copying, which is in conflict with the time invariance of that measure under 
Liouville dynamics. 

The relative entropy of the probabilities defined by Kullback and Leibler Dkl{T j i,V2) = 
J dxVi ln^i/T^) is constant in time when the probabilities obey the common Liouville 
dynamics [5,6]. This implies that the Kullback-Leibler distance, i.e., the relative entropy, 
remains the same as that given by the initial probabilities with which the original system 
is prepared. The Csiszar f-divergence [7] measures the distance between two probability 
distributions V\ and Vi by the form T = J dxVif^Pz/Vi), where / is a convex function 
satisfying /(l) = 0. Then, the Kullback-Leibler distance is found to be a special case of the 
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f-divergence (/ = -ln[V 2 /Vi], I 8 ])- 

As an extended proof of the classical no-cloning theorem along the line in [5], Plastino 
and Daffertshofer [9] have showed that the same contradictory situation upon copying a 
source system arises when employing a wider class of distance measures, i.e., the Csiszar 
f-divergence. Under the Liouville dynamics, the conservation of the f-divergence is incom- 
patible with a change in the distance during the copying process, where the following three 
different reasonable premises for the input and output distributions are made: 

(i) For the initial distributions, the state of a system is expressed by a product of distribu- 
tion functions of three components as Vj = J>^'p( a ^'p( r \ This reads as follows: the state 
of the original system (o) is copied to the target system (replica) (r) [10] with the help 
of an ancillary copying machine (a). The index j specifies the different initial conditions. 
The coordinates comprise the support of each distribution, i.e., V^(x^). This implies 
that the original system to be copied is prepared in such a way that its state is statistically 
independent of the machine and the replica. 

(ii) The successful cloning in (i) yields the final distribution as Qj = Q^Vj°\ During 
the copying process, the machine interacts with both the original system and its replica, and 
it functions as a reader and a writer. As a result, it changes its states from to Q( a \ 

(iii) The marginalization for the final distribution Qj with respect to the coordinates of the 
machine is used: / dx^Qj = V { ° ] {x^)V { ° ] {x^ r) ). 

With this background, we ask in this Letter: is the proof of the existence of the classical 
no-cloning theorem specific for the choice of the relative entropy? There may be other 
relative distance measures that also display a conflict under Liouville dynamics. To answer 
this question we investigate the case of non f-divergence-type distances with which, if valid, 
we can provide robustness to the delineation of the classical counterpart of the quantum 
process. We show that the combination of the f-divergence-type distance for measuring 
probabilities and the Liouville dynamics is a strong requirement for establishing the classical 
counterpart. We exemplify it by considering non f-divergence distances under a Liouville 
type equation incorporating a power law nonlinearity. In the next section, we first introduce 
a simple class of distance measures that do not belong to the f-divergence class. In section 

III, we show the discrepancy that leads to the equivalent conclusion given in [5,9]. In section 

IV, we consider nonlincarities in the Liouville equation for two types of distance measures, 
before closing with a summary and concluding remarks. 

A distance of non-Csiszar f-divergence type. — We consider the continuous Cher- 
noff a-distance [11] between two distinct probability distributions under phase space dynam- 
ics x t , 



C Q (Pi(f t ),P 2 (x t )) = -log 



J dx t Vf{x t )Vl- a {x t ) 



(1) 



We exclude the situation V" V\~ a = 0, that is, the two distributions do not have a common 
support, where the distance is not well defined. From an inequality a a b 1 ~ a > min{a, 6} 
for (0 < a < 1), the positivity C a (Vi,V 2 ) > is assured allowing for using it as a dis- 
tance. Similarly, it is asymmetric with respect to V\ and V2 for general value of a. The 
distance is zero if and only if V\ = V 2 holds for a fixed value of a. When V\ 7^ 7^2, 
then Cq{Vi,V2) = Ci(Vi,V 2 ) — holds. The Kullback-Leibler distance is generated from 
the differential coefficient, e.g., D KL (Vi,V2) = -[dC a (Vi,V2) /da] a =i and D KL (T , 2,T , \) = 
[dC a (Pi,P2)/da] a= o. A particular choice a = 1/2 provides the Bhattacharyya distance [12], 
C 1/2 {VuV2) = - log [J dxtVWP^- 

Proof of the no-cloning theorem. — Our probability distribution V(xt) is assumed 
to evolve according to the Liouville equation dtP + V • (vP) = 0, where v = dxt/dt denotes 
the velocity associated with the generalized coordinates in the phase space. Then, the 
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derivative of C a yields 
dC, 



dt 



= -\JdS t (aVr'Vt^ + (1 - a)V?V 2 - a ^ 



= ^ J dx t 



a 



l-a 



Vim) + (1 - a) (Jj^j V{vV 2 



7? 1 + (l-a)V(^) P_> 



dx t 



(2) 



where we put z = J dxtPf V 2 ~ a . In the second line, we substituted dtPj = —V • (vPj), 
(j = 1,2). The first term in the last line vanishes due to the Gauss theorem. Further, the 
calculations show that the second term vanishes. Therefore, the time derivative vanishes, 
i.e. dCa/dt = 0. Following the line of argument addressed in [5,9,13], we discuss the 
particular forms of states before and after the cloning process. We need to consider a 
tri-partite system whose distribution function is assumed to be factorized into component 
distributions as mentioned in (i). By the definition provided in Eq.(l), for the two factorized 
distributions for systems a and b, the total distance between two probability distributions V\ 
and V 2 satisfying Eq.(l) is equivalent to the sum of the distances between the two systems 
a and b, 



C a (Vx (o)Pi (6) , V 2 (a)V 2 (b)) = C a {P 1 (a) , V 2 (a)) + C a {V x (b) , V 2 (b)) . 



(3) 



Therefore, for the initial ensemble Vj = - pj°* l 'p( a )'p( r ) (j = 1,2), the distance is attributed 
to the distance between the original states 



c a {v u v 2 ) = c a (v[ o) A°\ 



(4) 



because we have C a (T (a \ V {a) ) = and C a {V {r \ V^) = 0. The two premises (ii) and (iii) 
mentioned in Introduction provide the contradiction in question. Indeed, the final distance 
in the form, Q, = Q (a) -pj o) , becomes 



C a {Q u Q 2 ) = 2C a {v[ o) ,vi o) ) + C a (Q\ a >, Q%>) 



■.(a) n (a)-, 



(5) 



This leads to a negative value of C a if we require that C a (Qi, Q 2 ) — C a (Vi,V 2 ) should 
hold. This, however, is in conflict with the positivity property of C a . 
Further, when we evaluate the C a (Qi, Q 2 ) based on the premise (iii), we have 

CaiQuQt) > -log^d^) jdtipQi Jda!i a) Qi 

= -2io g (y ' dx { :\v[ o) r{vi o) f-^ 



= 2C a (J> x ,V 2 ) 



(6) 



where the first inequality follows from the Holder's inequality; the second, from (iii); and 
the last line, from Eq.(4). The distance between two different states after copying becomes 
larger than or equal to the twice of the initial one. Hence, in either case, we could show 
the contradiction with Eq.(2) as long as V\ ^ V 2 holds for all values of a. In the present 
setting, the Liouville dynamics can not allow for classically copying an unknown original 
state to another system without destroying the original system. 
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Consideration of distance form in terms of a Liouville-like nonlinear evolu- 
tion equation. — The Liouville equation must be linear and governs the evolution of a 
probability density that describes a statistical ensemble (in the sense of Gibbs) of dynamical 
systems, all evolving according to the same equations of motion. Such a time-dependent 
ensemble probability density always evolves according to a linear equation (e.g. [14]). Start- 
ing from the Liouville theorem, which states the conservation of the probability distribution 
function along the orbit in the phase space, the derived Liouville equation cannot have pow- 
ered probabilities such as P q , where q may be relevant to deviation from linearity. However, 
it has been pointed out that the linearity of the Liouville equation does not preserve statis- 
tical independence when the linear combinations of probability distributions Q = cj Qj 
are applied to the associated final states with the marginality property [5] . In order to gain 
more insights into the relation between the dynamics of the classical statistical ensemble 
distribution in copying process and the distance we measure, we here consider two cases 
whose dynamics evolve according to a Liouville-like nonlinear evolution equation. The same 
reasoning for the classical no-cloning consideration should provide constraints for the func- 
tional form of the distances and the parameter of the dynamics upon copying. 

There can be two immediate generalizations of the Liouville equation that break its 
linearity, d t V q + V(vP) = and d t V + S7(vP q ) = 0. These forms are reminiscent of the 
nonlinear Fokker-Planck equations, where the probability V is powered for terms in the 
equations [15]. Since the former is equivalent to the latter by a transform V q — > V, we 
use the latter in the following consideration. It would be worth mentioning the differences 
between the standard Liouville equation and its nonlinear forms. The Liouville-like nonlin- 
ear equations may govern the behavior of a time-dependent density describing a set of real, 
interacting systems instead of an evolution of a probability density associated with a statis- 
tical ensemble. Indeed, the nonlinearity in the evolution equation for the density constitutes 
an effective description of the interaction between the systems and has relevant applications 
(e.g. [16]). In this sense, the nonlinear Liouville-like equations in our consideration are akin, 
for instance, to the Vlasov-Poisson equation (e.g. [17] and references therein). In spite of this 
difference, we believe that the use of the nonlinear Liouville-like equation in the followings 
serves as a supporting tool for our present discussion. 

First, we consider a distance G\ defined by the following form 



where h is a function of the ratio of two different distributions r\ = VijV\. When a = 1, 
Gi reduces to the Csiszar's f-divergence type, which has already been considered in [9]. We 
are interested in how the exponent a governs the form h when a ^ 1 and a^O. Under the 
evolution, which can be described by a nonlinear generalization of the ordinary Liouville 
equation, 



where q is not equal to unity, we consider the constraint on the functional form of h that 
makes the distance G\ time independent. We have 



where h' means the derivative with respect to n. Putting the time derivative of Vj as 
—V(vPj) from Eq.(8) and from the Gauss theorem, we obtain the following expression 




(7) 



d t V + V{vP q ) = 0, 



(8) 




(9) 




p-4 



Title 



- J v [PfV [V?-\ah - V h')] + VlViV^h')] dx t . (10) 
For the rhs. of Eq.(9), we find 

Jv[(a- 1)V?- 2 {aVfh - ti( V Vl - )} VVi 
+ V?- 1 {{a - l)Vfti - (riVf -V2)h"}Vr]]dx t . (11) 
Because the coefficients of VPi and Vrj have to vanish, we obtain 



aVfh = (rjVl - Vl)bl 
{a-\)V\hl = (f ] Vf-V2)h". 



(12) 



Simultaneous satisfaction of both relations in Eq.(12) leads to the ordinary differential equa- 
tion b! = ch^i with the integral constant c, which is of the Bernoulli type and can be solved 
by the change of variables z — h 1 ^ . Together with the condition h(l) = 0, we have 



l-u 

1-a 



h ^ = [—) (v-iy- a - as) 

Note that the nonlinear parameter of the dynamics q does not enter into the functional 
expression. Conversely, Eq.(13) must satisfy Eq.(12). Substituting Eq.(13), e.g., into the 
first equation of Eq.(12), we have for any i] 



a rj — rf 



1 — a rj — 1 1 



(14) 



therefore, we can conclude q = and a = 1/2. This case, however, is trivial, indicating that 
the combination of Eq.(7) and Eq.(8) may be applicable only to the case q = 1. 

Let us now consider the case a = 1 in terms of another possible distance measure. 
We see that dG/dt = can be achieved robustly under the ordinary Liouville equation 
dtP + V(vP) = by investigating the following combination: 

G 2 = j dxtV^m, (15) 

where (3 is kind of a weighting parameter for the second probability distribution. The similar 
calculation for this provides the result 



dG 2 



/dt = Jv [f3((3 - l)V q 2 +f3 - 2 (VV 2 )f + (pVt - 1 - »?7>?)/"Vrj] dx t . (16) 



The necessary conditions of dG 2 jdt = are (3 = or (3 = 1 from the first term. However, 
f3 = implies that G 2 does not contain V 2 , which should be excluded from the present 
consideration. From the vanishing second term, we obtain a relation j3 — (V 2 /Vi) 1 ^ q . 
Therefore, as long as V\ 7^ V 2 holds, [3 = 1 yields q = 1. 

Summary and concluding remarks. — We have shown that the universal cloning 
machines are incompatible with the conservation of information distance that is different 
from the Csiszar's f-divergence under the Liouville dynamics. The Chernoff a-divergence 
provides an example that makes the cloning under that dynamics infeasible. That is 
dCa/dt — holds in spite of C a (Qi, Q 2 ) ^ C a (Vi,V 2 ) before and after the copying process. 
One question still needs to be answered: to what extent the present argument is valid? We 
answer it as follows. Since a general form Q = J dx t Vif(V 2 /Vi) satisfies dG/dt = under 
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the usual Liouville dynamics [9], we can see that the discussion developed here is also true 
for the more generic distance measure G = G{Q) because it immedeately provides 

dG dQdG n 

lE-mog- - (17) 

From this perspective, the Chernoff a-distance we employed is only a special case (when the 
function G is logarithmic) that makes the classical counterpart of the quantum no-cloning 
theorem hold. 
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